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Abstract 

We investigate the attractor mechanism for spherically symmetric extremal 
black holes in a theory of general R 2 gravity in 4-dimensions, coupled to gauge 
fields and moduli fields. For the general R 2 theory, we look for solutions which 
are analytic near the horizon, show that they exist and enjoy the attractor 
behavior. The attractor point is determined by extremization of an effective 
potential at the horizon. This analysis includes the backreaction and sup- 
ports the validity of non-supersymmetric attractors in the presence of higher 
derivative interactions. To include a wider class of solutions, we continue 
our analysis for the specific case of a Gauss-Bonnet theory which is non- 
topological, due to the coupling of Gauss-Bonnet terms to the moduli fields. 
We find that the regularity of moduli fields at the horizon is sufficient for 
attractor behavior. For the non-analytic sector, this regularity condition in 
turns implies the minimality of the effective potential at the attractor point. 
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1 Introduction 



Supersymmetric BPS black holes in string theory are known to exhibit an attractor 
mechanism, whereby the values of scalar fields at the horizon are determined only in 
terms of the charges carried by the black hole and are independent of the asymptotic 
values of these scalar fields [1-3]. As a result, the black hole solution at the horizon 
and the resulting entropy turn out to be determined completely in terms of the 
conserved charges associated with the gauge fields. Following this, the attractor 
mechanism has been studied extensively in the context of the supergravity and 
string theory. 

It was further discussed in [4,5], that the attractor mechanism can also work for 
non- supersymmetric cases. Recently there has been a surge of interest in studying 
the attractor mechanism without the use of supersymmetry. In particular, impor- 
tant developments are taking place in the context of non-supersymmetric attractor 
mechanism, in the presence of higher curvature interactions in theories with and 
without supersymmetry. 

The role of non-supersymmetric attractors was further clarified in [6]. By using 
perturbative methods and numerical analysis it was shown that the attractor mecha- 
nism can work for non-supersymmetric extremal black holes. The authors considered 
theories with gravity, gauge fields and scalars in four and higher dimensions which 
are asymptotically flat or Anti-De Sitter. It is worth mentioning that the scalar 
fields do not contain any potential term in the action, so as to allow them to behave 
as moduli at infinity. However, the coupling of these scalars with the gauge fields 
acts like an effective potential for the scalars. This is with the foresight to fix the 
scalar field values in terms of the conserved charges associated with the gauge fields. 
Varying the effective potential, the values of scalar fields are determined in terms of 
the charges carried by the black hole. Several examples were presented in [6] and 
also in the context of string theory [7]. Recently in [8], a c- function was introduced 
which monotonically decreases from infinity to the horizon and coincides with the 
area of the horizon at the horizon. Further developments are studied in [9-11]. 

Following these developments, it is important to understand whether the attrac- 
tor mechanism can work in the absence of supersymmetry when there are higher 
derivative terms in the action. In [12], using the near horizon geometry of extremal 
black hole to be AdS-z x S 2 , it was explicitly shown that extremizing the 'entropy 
function' with respect to the scalar fields, leads to a generalized attractor mechanism 
(see also [13]). The Legendre transform of this function with respect to the elec- 
tric charges, evaluated at the extermum, is proportional to the entropy of the black 
hole. This method has been used to compute corrections to the entropy of different 
black holes due to the higher corrections to the effective action [14-17]. In fact, 
the addition of a Gauss-Bonnet term reproduces the entropy associated with four 
charge black holes in the Heterotic string theory, which were earlier computed mak- 
ing heavy use of supersymmetry. All these results also agree with the microscopic 
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counting of entropy [18-28]. However the existence of a solution is a presumption 
in [12], though it has the advantage that the action is very general. 

Recently it has been shown that the Legendre transformation of black hole en- 
tropy with respect to the electric charges is related to the generalized prepotential, 
and this led to the new conjectured relation between the black hole entropy and 
topological string partition function [29-32]. Applying the results of these black 
holes to special case of black holes in Heterotic string theory with purely electric 
charges, one finds agreement between black hole entropy and the degeneracy of el- 
ementary string states [33-41] even for small black holes with vanishing classical 
horizon area [42-44]. 

Following the approach in [6], we show that the attractor mechanism works for 
non-supersymmetric theories even after the inclusion of higher derivative R 2 terms, 
and in a more specific case, the Gauss-Bonnet terms. In four dimensions, since 
the Gauss-Bonnet term is a total derivative, its coupling with dilaton field plays an 
important role. Without this coupling, Gauss-Bonnet term dose not contribute to 
the equations of motion. Thus, after coupling to the scalar fields we find non-trivial 
equations of motions and solve all of them together. 

To solve the equations, we consider the background to be an asymptotically flat 
extremal black hole and look for scalar fields which are free at infinity and regular at 
the horizon. It will be shown that subject to these boundary conditions, all solutions 
of the scalar fields get fixed values at the horizon, independent of their arbitrary 
moduli values in the asymptotic region. This analysis supports the existence of 
attractor mechanism for black holes in higher derivative gravity. 

Since we are faced the situation of dealing with a set of coupled non-linear equa- 
tions involving the metric and many scalar fields, we solve them by a series solution 
method. In R 2 gravity we derive the solutions as a series with analytic behavior near 
the horizon. For completeness, we consider non-analytic, but regular solutions in 
the Gauss-Bonnet case. We also make use of a perturbative method in terms of the 
parameter a', for the case of Gauss-Bonnet theory and construct the perturbative 
solutions on top of the zeroth order solutions of the Einstein-Hilbert gravity given 
in [6]. In all these cases, we show that, as long as regularity and asymptotic flatness 
conditions are satisfied, the solutions enjoy the attractor mechanism. 

The organization of the paper is as follows. In the next section we give a brief re- 
view of the non-supersymmetric attractor mechanism in the Einstein-Hilbert theory. 
In section 3, we discuss the attractor mechanism near the horizon of an extremal 
black hole, for general higher derivative theories. Section 4 is devoted to R 2 gravity 
where we find an analytic series solution for the set of coupled equations. In section 
5, we investigate Gauss-Bonnet theory and give as large as possible class of solutions 
with given boundary conditions. We conclude in section 6. 
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2 Brief Review of Non-Supersymmetric Attrac- 
tor Mechanism 

Here we collect some important points regarding the non-supersymmetric attractor 
mechanism in four dimensional asymptotically flat space-time for later use. It is 
instructive to start with gravity theories coupled to U(l) gauge fields and scalar 
fields, dictated by following the bosonic action 5 : 

S = i J d'x^g (r - 2(<%) 2 - f ab (<p l )F^F b r - -^f^)F^F%€r^ 

9 (1) 

where F® u , a = 0,...N are gauge fields and <p l , i = l,...n are scalar fields. The 
scalars have no potential terms but determine the gauge coupling constant. We 
note that refers to the metric in the moduli space and this is different from space 
time metric, i.e., g^ v . The equations of motion derived from the action (JTJ are as 
follows: 

itV - Id^Mi = fMi) (2F° M F 6 / - \g^F a KX F b " x 

■ -MV^d^) = -MU)F*F h ^ - \*d i {u)nj%<r>° i 



d, U^gfaMF^ + \f ab F b p(J e^ = 0. (2) 

A spherically symmetric space-time metric in 3+1 dimensions can be taken to be of 
the form: 

ds 2 = -a{rfdt 2 + a{r)~ 2 dr 2 + b(r) 2 dtt 2 (3) 

On the other hand, the Bianchi identity and equations of motion of gauge fields can 
be solved by taking the gauge field strengths to be of the form: 

F a = f a \&) (Qeb - fbcQ c J ^dt A dr + Q a m sin OdO A dip, (4) 

where Q a m and Q ea are constants that determine the magnetic and electric charges 
carried by the gauge fields F a , and f ab is inverse of f a b- Using the equations of 
motion it is possible to derive the following second order equations for the unknown 
functions a(r), b(r) [6]: 

{a 2 {r)b 2 {r))" = 2, (5) 
~{dA)\ (6) 



d r (2a 2 b 2 d r ^) = (7) 
l+a 2 b /2 + ^~ = -^(K//(^)) + aV(^) 2 , (8) 



3 In the convention of [6], ^= factor is involved in the definition of e^ vpa . 
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with the effective potential given by, 



Vefffa) = f a \Qea ~ facQ^Qeb ~ f bd Q d J + fabQ^Q^ (9) 

It is useful to note that the equations of motion (jUHZj) given above, can be derived 
from the following one-dimensional action: 



S = I 



2,2\" o.2,," o.2,2/o ,n2 2V eff( 



(a 2 b 2 ) -2a 2 bb -2dW(d r <t>y 



b 2 



(10) 



Furthermore, eq. (jHJ) stands for the Hamiltonian constraint which must be imposed 
in addition. 

Now two sufficient conditions having the attractor behavior for the moduli fields 
can be stated as follows [6]. First, for fixed charges, V e ff as a function of the moduli, 
must have a critical point. Let us denote the critical values of the scalars as (j) 1 = </>q. 
Then we have, 

0iK//(&o) = o. (11) 
Second, the matrix of second derivatives of the potential at the critical point, 

M tj = ^didjVeffifa), (12) 

should have positive eigenvalues. Schematically one can write, 

M i:j > 0. (13) 

This condition guarantees the stability of the solution. Let us refer to as the 
mass matrix and the corresponding eigenvalues as masses (more correctly mass 2 
terms) for the fields, <j) [6]. 

Once the two conditions mentioned above are met, it was argued in [6] that 
the attractor mechanism typically works. There is an extremal Reissner-Nordstrom 
black hole solution in the theory, where the black hole carries the charges specified by 
the parameters, Q^, Q ea . Further, the moduli take critical values, <pQ at the horizon, 
which are independent of their values at infinity. In other words, although <j) is free 
at infinity as a moduli field, it is attracted to a fixed point, </>o at the horizon. 

Let us now consider a constant 0o solution, in which case, the resulting horizon 
radius can be shown to be, 

b 2 H = Vtffifa), (14) 

with the entropy given as: 

Sbh = \A = irb 2 H (15) 

In the next section, we study the behavior of solutions near the horizon and discuss 
the attractor conditions for more general theories. 
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3 Near horizon solution 



Consider, a general theory of gravity and gauge fields with higher derivatives inter- 
actions. Let us assume that the theory admits an extremal Reissner-Nordstrom(RN) 
black hole solution as in (J3J). Let us consider the general action as, 

S = J d 4 xy/g C{g^, v , F^, //), (16) 

where //s are some couplings in the theory. To include scalar fields, we replace 
these couplings by some functions of scalars (but not their derivatives) as fi(4>), and 
add a kinetic term for the scalars in the action. Besides the modifications to the 
equations of motion of metric and gauge fields, one also has the equations of motion 
of scalar fields which take the form, 

r% (aV3^) = (17) 

where W e ff comes from the effective Lagrangian density, after inserting the ansatz 
for the metric and gauge field solutions as background fields in the action. 

The near horizon behavior of this equation is important. To a first order approx- 
imation, we neglect the backreaction and assume the existence of a black hole with 
double degenerate horizon. In this case, one expects to have the following forms for 
a(r) and b(r) near the horizon: 

a(r)~(l-^), 

6(r)~r. (18) 

Now, for a regular solution of 0(r) near the horizon, the left hand side of ()17|) 
vanishes, which in turn implies that W e /f should have a critical point as r approaches 
Th- Let us denote this critical point by <frio, such that: 

■*fi*=°- (19) 

Notice that in contrast to the previous section where V e ff was r-independent, 
here W e ff generically depends on r, but its extremum will be meaningful only at 
r — th- Thus, we use the subscript H which indicates the value of W e ff at the 
horizon. 

In the other limit at large r, if we assume an asymptotically flat background, then 
the scalar field should appear as a moduli field with arbitrary asymptotic values. 
In this way, one can obtain the attractor mechanism, which means that, for any 
arbitrary values of moduli fields in the asymptotic region, the scalar fields always 
approach fixed critical values (f)^ at the horizon. 
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In [6], this behavior was studied by the perturbative methods described above. 
That is, by setting the asymptotic values of the scalar fields equal to their critical 
values and then examining what happens when the scalars take values at asymptotic 
infinity which are somewhat different from their attractor values, (pi = (piQ. 

To this end, one starts with a study of the scalar field equation to first order in 
perturbation theory, in the RN geometry without including backreaction. Let (pi be 
a eigen mode of the matrix of second derivatives, = didjW e ff. Then, denoting 
S(pi = (pi — (pio, neglecting the gravitational backreaction and working to a first order 
in S(pi, we find that eq.([17|) takes the form, 

d r ((r - r H ) 2 d r 5(j)i) = ■i r w 2i 5(j) i , (20) 

where W2i is the relevant eigenvalue of Wij. In the vicinity of horizon, we have 
replaced b(r) by r#. Asymptotically, as r -> oo, the effects of the gauge fields 
and higher derivatives interactions die away and eq.([2L)jl reduces to the equation of 
motion for a free field in flat space. This has two expected solutions, 5(pi ~ const. 
and 5(pi ~ 1/r, both of which are well behaved. One can now notice that the second 
order differential equation is regular at all points in between the horizon and infinity. 
So once we choose the non-singular solution in the vicinity of the horizon it can be 
continued to infinity without blowing up. 

Now one can easily find the solutions of the linearized equations (J2*U|) . to be [6]: 

Mr)=c(l-H) ) * (21) 



/• 

where, 




(22) 

With a plus sign, this represents a regular solution at the horizon, provided W2i is 
non-negative. Further, in eqn.(21) C is an arbitrary constant, which denotes the 
asymptotic value of the scalar field. 

Next, one includes the gravitational backreaction. The first order perturbations 
in the scalars source a second order change in the metric. The resulting equations 
for metric perturbations are regular between horizon and infinity. The analysis near 
the horizon and at the infinity shows that a double- zero horizon black hole solution 
continues to exist which is asymptotically flat after including the perturbations. In 
case of [6] the analysis was extended to all orders in perturbation theory analytically 
and it was found that the attractor mechanism works to all orders without any extra 
conditions. 

To conclude, the key feature that leads to the attractor behavior is the fact 
that solutions to the linearized equation for 5<p are well behaved as r — »• oo, and 
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under some conditions there are regular vanishing solutions near the horizon. If one 
of these features fails then the attractor mechanism may not work. For example, 
adding a mass term to the scalars, results in one of the two solutions at infinity 
diverging. Now, as discussed in [6], it is typically not possible to match the well 
behaved solution near the horizon to the well behaved one at infinity and this makes 
it impossible to turn on the dilaton perturbation in a non-singular fashion. In 
the following sections we consider some interesting examples to study the attractor 
mechanism. 

Before going through the examples, we wish to make a comment about the en- 
tropy of the black hole. Indeed, following the Wald's formula or using the prescrip- 
tion for the entropy functional given in [12], one finds the entropy to be proportional 
to W H (cj) ) (see [12,14,16,17]). 



4 R 2 - Gravity 

In this section, we consider a general R 2 action coupled to the moduli field and add 
it to the Lagrangian (£[]): 6 

C 2 = G(0) ( V R 2 + PR^BT + SR^R^) . (23) 

The most general metric with spherical symmetry is of the form: 

ds 2 = -a(r) 2 dt 2 + . J . dr 2 + b(r) 2 d£%. (24) 

Although it is possible to set c(r) = 1 by a redefinition of the radial coordinate, we 
keep c(r) free to derive the complete set of equations, including the Hamiltonian 
constraint, from the one dimensional action. 

The gauge field background is same as in eqn. (j3J). Inserting the metric (J2"4^ in 
the action, we find the following one dimensional action 7 : 



S = [dr --- (a 2 b 2 )" - ca 2 bb" - _ 2 ba 2 b'c' - b 2 ac'a' 

J [c 2 K ' b 2 c 

+^G(4>) [2r] (4ac 2 a'bb' + cb 2 a'c'a + a 2 c 2 b' 2 + c 2 b 2 aa" + 2bca 2 b'c' - 1 

+cW + 2bc 2 a 2 b") 2 + (3(1- 2a W + c 4 oVa" 2 + 3c 4 6 2 6"V - 2bca 2 b'c' 
+ 8cWa 3 a'b' + 8c 3 b 2 b' 2 c'a 3 a' + 2bc A a%"b' 2 + 4a 3 c 4 aW 3 + 3c 2 b 2 b' 2 c' 2 a 4 
+2c 4 b 3 a 3 a"b" + 2c 3 b 3 a 3 a"b'c' + 2c 3 b 3 a'c'a 3 b" + 2c 2 b 3 a'c' 2 a 3 b' + 2c A b 3 a' 2 b"a 2 



6 For simplicity we consider a single moduli field. The generalization to multi-moduli fields is 
straightforward. 

7 This is the same as the case of action (I1UH where beside the V e ff term, all the other terms (being 
gauge-background independent) are simply derived by inserting the metric ansatz into the original 
action. V e ff term can be reached by inserting the gauge field solution with a sign modification. 
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+ 6c 3 b 2 b"a 4 b'c' + 2c 3 b 4 a 2 a"a'c' + 2c 4 b 4 aa"a' 2 + Ac 4 b 3 a 2 a" a'b' + c 2 b 4 a' 2 c ,2 a 2 
+2c 3 b 4 a' 3 c'a + 6cWV 'a 2 b' + Ac 4 b 3 a' 3 ab' + 8c 4 b 2 a 2 a' 2 b' 2 + c 4 b 4 a' 4 + a 4 c 4 b' 4 
-2bc 2 a 2 b" - Aac 2 a'bb' + 2a 4 c 3 b' 3 bc') + 25 (l + 2c 3 b 4 a 2 a"a'c' + 2c WV 
+c 4 b 4 a 2 a" 2 + Ac 3 b 2 b"a 4 b'c' + c 2 b 4 a' 2 c' 2 a 2 + Ac 4 b 2 b"a 3 a'b' - 2a W + c 4 6 4 a' 4 
+ 2cWc'V + Ac 3 b 2 b' 2 c'a 3 a' + 2c 4 b 4 aa"a' 2 + 2c 3 b 4 a' 3 c'a + Ac 4 b 2 a 2 a f2 b' 2 
+a 4 c 4 b' 4 )}] (25) 

One derives the equations of motion by varying the above action with respect to 
a, b, c and <fi. Then we can put c(r) = 1, and the equation for c turns out to be 
the Hamiltonian constraint. We call these equations as EqA, EqB, EqC and Eq§, 
respectively 8 . Though the total effective potential W e ff introduced in the previous 
section has a complicated form in this theory, it is not too hard to find its critical 
point at the horizon as we show below. 

It is well known that these equations admit AdS^ x S 2 as a solution with con- 
stant moduli. However we want to address the attractor behavior in solutions 
which are asymptotically flat. In view of the fact that the four equations governing 
(a(r), 6(r), 4>{r)) are a set of highly complicated coupled differential equations of or- 
der four. To solve these equations, we follow the Frobenius method. As a variable 
of expansion we define x = 1 — — , ranging from to 1 to cover r > r# completely. 

In fact, requiring that the solution: 

1. be extremal: meaning that a(r) = (r — r^) 2 i/(r), with H(r) being analytic at 
the horizon, r = r#. 

2. be asymptotically flat: meaning that (a(r), 6(r), 0(r)) tend to (l,r, 0^) at 
asymptotic infinity. 

3. be regular at the horizon, 

the most general Frobenius expansions of a(r), 6(r) and <j)(r) take the form: 

oo 

a(r) = {x + x Xl ^a n x n ), (26) 

n=2 

oo 

6(r) = (l + x A ^ &„*"), (27) 

(1-*) ^ 

oo 

0(r) = (0 o + x A3 Yl <f>nx n ), (28) 



n=l 



with A,- > 0. 



For concreteness, from now on we set both the functions G(<f)), which is the 
coupling of moduli field to R 2 , and V(<f>), which is the U (l)-background contribution 

8 Here we avoid writing these equations which are very lengthy. We have used the maple package 
to derive them. 
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to the effective potential, to be of standard dilatonic forms, as inspired by string 
theory. More precisely, denoting the attractor value of moduli with </>o we set, 9 

G(<j>) = g e a ^~ M ; V(<j>) = v e ^~M. (29) 

We mention that although, in comparison to (J9~|). V(<p) of (J2*9")l is of pure magnetic 
(electric) type, the case given in (|2T)|) is rich enough for our study since the new 
function G(0), mathematically, compensates the role played by the lacking electric 
(magnetic) term. 

Zeroth order results 

For the case at hand, it is consistent to set A« = 0. EqA is automatically satisfied 
and either of EqB or EqC implies, 

v = r 2 H , (30) 

while Eq§ gives: 

<™> = ■ (31 » 

Notice that equations (j3*Uj) and ()31|1. together, determine both of the attractor value 
of the moduli filed and the horizon radius in terms of the parameters of the action. 
In fact both the above results are meaningful to us. Due to ()3()jl the Bekenstein- 
Hawking entropy of the solution is given by the value of the parameter vq, up to 
a numerical prefactor. The information encoded in (}3*T]) becomes clearer when we 
note that this relation, together with (|2~9"j). implies: 

§ -2(/3 + 25)§ = 0, (32) 
d<p acp 

which is exactly ()19j) with, 

W H {(t>) = V{(j>) - 2{f3 + 28)G((j)). (33) 

This in fact fixes 0o & t its extremum point. From (J28)) . 0o — 0( r n) and so the value 
of the moduli field is fixed at the horizon, regardless of any other information. Thus 
to complete the proof of the attractor behavior, we should be able to show that the 
four sets of equations of motion, denoting a coupled system of differential equations, 
admit the expansions (|26|) . ()27)1 and (J2EJ)- Furthermore, one should see that there 
are solutions to all orders in the x-expansion with arbitrary asymptotic values at 
infinity, while the value at the horizon is fixed to be 0o- The existence of a complete 
set of solutions with desired boundary conditions (considering the fact that we have 
coupled non-linear differential equations) by itself is not trivial. Moreover, it is easy 
to show that, in our theory, there is no asymptotically flat solution with everywhere 
constant moduli. 

9 Notice that with this parametrization, G(4>q) = go and V{4>o) = v o- 
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Figure 1: 0(r) vs. log(r), where numerical coefficients are r] = 1, /3 = 2,8 = 3, fi = 5 
and go = 4. Different curves represent different asymptotic values for K = 0^. The 
attractor point is 0o = 1 at the horizon, rjj = 1. 

Higher order results 

We solve the system of equations (EqA, EqB, Eq$, EqC) order by order in the x- 
expansion. To first order, we find that one variables, say <fii, can not be fixed by the 
equations. Let us denote the value of 0i as K. We thus find a 2 and b\ as a function 
of K. One can check that at any order n > 2, one can substitute the resulting values 
of (a m , 6 m _i, c m _x), for all m < n from the previous orders. Then (EqB, Eq§, EqC) 
of the current order together with EqA of order (n — 1), consistently give, 

b n = b n (K,g ) ; a n+1 = a n+1 (K , g ) ; (j) n = (j) n (K, g Q ) . (34) 

as polynomials of order n in terms of K. 

It is worth noting that K remains a free parameter to all orders in the x- 
expansion. From (J2TI)) . ([2~7j) and the asymptotic values of (a{r),b{r),cf){r)) 

are given by a sum of all the coefficients in the x-expansion of the corresponding 
function. Owing to the result ()34|1 we observe that (a^, b^, 4>oo) are f ree to take 
different values, given different choices for K 10 . 

The arbitrary value of <fi at infinity is = 0oo, while its value at the horizon is 
fixed to be <J)q. This is the very attractor mechanism which we were looking for in 

10 We do not address the convergence of the series in detail, but it would be the case for small 
enough values for \K\. 
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this model. Figure 1 shows <p vs - r m logarithmic scale with different asymptotic 
values K. 



5 Black Holes in Gauss-Bonnet Gravity 

Since the Gauss-Bonnet case is a special sector of general R 2 gravity, it is obvious 
that any solution of the form considered in the previous section, enjoys the attractor 
mechanism. Now we relax the analyticity condition and, instead, only require the 
solutions to be regular at the horizon. This relaxation, enables us to cover a class of 
solutions connected to G = case discussed in [6], which are not analytic at x — 0. 

Black holes in higher order gravity are being actively pursued. It is important 
to note that string theory predicts the Einstein-Hilbert action to be modified by 
higher order curvature corrections. The leading correction is quadratic in curvature. 
Typically, only certain special corrections are considered. At quadratic order, the. 
combination which is considered is the Gauss-Bonnet term: 

TZ jS XaR lSXa - IRjsR 15 + R 2 (35) 

We are interested in an Einstein-Gauss-Bonnet theory coupled to the gauge fields 
and scalars. Since the Gauss-Bonnet term in four dimensions is a topological term, 
the coupling to the scalar fields is crucial. Without this coupling, the Gauss-Bonnet 
term does not contribute to the equations of motion. In that case, we are left with 
the RN solution reviewed in section 2. 



5.1 Equations of motion and attractors 

Consider a gravity theory coupled to U(l) gauge fields and scalars in the presence 
of higher order corrections of the Gauss-Bonnet type. The system is dictated by 
following bosonic action: 

S = 1 J d'x^g^R - 2(«%) 2 - f ab ^)F%F b ^ - ^f ab ^)F^F b pa e^ 

+aG(<j>) (R^xaR 15 ^ - 4R jS RP< s + R 2 )^j, (36) 

where a is the Gauss-Bonnet coefficient with dimensions of (length) 2 and is positive 
in Heterotic string theory. The function G(4>) is arbitrary and signifies a dilaton-like 
coupling. 

Let us also assume that all quantities are functions of r and take the follow- 
ing ansatz for the metric and gauge fields of a Reissner-Nordstrom- Gauss-Bonnet 
(RNGB) black hole: 

ds 2 = -a(r) 2 dt 2 + a(r)- 2 dr 2 + b(r) 2 dn 2 , (37) 
F a = f ab (<p l )(Qe b -f b cQ c J^dtAdr + Q a m sznedeAd<P, (38) 
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Plugging 11 the above ansatz in the action one finds the following one-dimensional 
action: 



5=4 j 



H 2 



,2r2V o„2ll" o„2l£/"fl -N2 2V effWi) 



2-(a 2 b 2 ) -2a 2 bb -2a 2 b 2 {d r <t>y 



b 2 

+ 8aG((f))d r (aa'(-l + a 2 b' 2 ))] . (39) 
The Hamiltonian constraint, which must be imposed in addition takes the form: 

_ i + a 2 fe '2 + + * Veff (fc) - a 2 b 2 ^ 2 = —AaG' (aa' - 3aV&' 2 ) , (40) 

where a primes indicates derivation with respect to r. 
The equations of motion following from ([39)1 are 12 : 

(W = ss^ + ^ K (i- v))', (4i) 

(a 2 6)"-a /2 6-2aa'6 / + aV 2 + ^ = 4a ((a 3 aW)' + a 3 aW) , (42) 

bb" + b 2 (j)' 2 = a(-2G" + 2a 2 {G'b' 2 )'y (43) 

To solve the above equations, we take two different approaches: re-expansion and 
a-expansion as explained in the following subsections. 

5.2 x-Expansion series solutions 

By ^-expansion, we mean a set of solutions as a Frobenius series given in (|26| ) -(|28j ) . 
We take a common A» = A for all the solutions and write the series as: 

<f>{r) = 0o + Kx x + ■■■ , 

a(r) = x + aiX A+1 + • • • , 

b(r) = r(l + 6ix A H ). 

Let us also consider Taylor series expansions for V(4>) and G(<p) as follows, 

V{4>) = w o + ^i(0-0o) + ^ 2 (0-0o) 2 + --" , (44) 
G(4>) = go + gi(4>-^o) + ^92{<p-<Po) 2 + ■■■ ■ (45) 



11 See footnote 7. 

12 It is easily checked that the RHS of the first equation, if evaluated at the horizon, leads to the 
effective potential W e f / = V + 4 a G. This concides with the general effective potential (I33|) for 
the G.B. parameters (/3, 5) = (—4, 1), if we re-absorb a into the function G. 
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Since integer values of A are included in the results of section 4, we look for non- 
integer values A > 0, which give us a regular (though non- analytic) set of solutions 
at the horizon. 

By a careful investigation near the horizon, for the lowest power of x which is 
x A , one can solve the set of equations together and find the non-trivial solutions as: 



with <7i = ai = b± — 0. However, i> 2 , gi and K are undetermined to this order. 
The value of A in (J46|) is precisely what we already found in eqn. (}22|) with Wh = 
V{(j>) + 4G(0). The second equation in (|4Tj) is the extremum condition for Wh 
which gives the attractor value 0o at the horizon. Notice that we are faced with an 
extra condition gi = 0, which indicates that G, V and Wh are at their extremum 
at the horizon, simultaneously. Such a case would be the only situation where a 
non-integer A can be found. Otherwise we have to choose A = for gi ^ 0. 

The regularity condition for indicates that A should be non-negative and it in 
turn gives Vi + Ag 2 > 0. This again means that Wh is minimum at its extremum 
point </>o- 

Higher order terms can be derived in a similar fashion. The important point 
is that, due to the non-linear nature of equations, they are a mixture of different 
powers of A, like x nX as well as x nX+m . To order these powers, we assume < A < 1. 
Then the next leading term would be x 2X . For higher order terms, since A is already 
known from the first order result ([46)1 . we can distinguish the next order from x 3X 
and x x+1 . For small enough A, it shows that we are generating a power series, x nX 
as argued in [6]. 

Notice that in contrast to the analysis of section 3, here we considered all the 
equations simultaneously. This first means that, in principal, we are taking the 
backreaction into account. Secondly, since we are dealing with a higher derivative 
theory, besides the so-called Klein-Gordon equation for field. Other equations also 
involve the second derivative of <fi and are important in the dynamics of 0. So, they 
should be investigated as well. 

5.3 a-Expansion 

Motivated by low energy effective actions of string theory, it is be reasonable to 
consider a as a small parameter (proportional to a') and try to solve the set of 
equations perturbatively in the a parameter. Obviously, to zeroth order, we start 
with the results in [6]. Let us consider 



A 




(46) 




(47) 



0(r) 



(o) + a0 (1) (r) + a 2 ( %) + ••• , 
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a(r) = (l - — ) + cra (1) (r) + a V 2) (r) + 

6(r) = r(l + ab {l) (r) + a 2 b {2) (r) + •••), 
together with Taylor expansions (}4*4*)) and (|45|). 



Take <pP> = 0o, the constant solution at zeroth order. Then we find, 



() = ( 1 -— ) + p 4(r), (4* 



r 



where P4 is a function of r which behaves like a polynomial x of order 4 near the 
horizon. Moreover a^\r) and ^(r) can be set to zero consistently. It is worth 
mentioning that the attractor equation (|47|) should remain valid to all orders in a, 
so that when expanding in powers of a: 

= vi+ (y 2 <j>$ + 4^) a + i (^g } + t^ff) 2 + 8^0^) a 2 + 0(a 3 ), 

= 0. (49) 
Solving at each order in a, one finds: 

«1 = 0, 



t(2) 



w 2 

-3 (2gi92V2 - g 2 v 3 ) 



where the first equation is the extremum condition to zeroth order in a and the 
other two equations give the correct boundary value of the first and second order 
fields at the horizon, respectively. 

At the second order, we one can find non-zero (r) , (r) and (p^ (r) , in terms 
of some integration constants. These constants are subject to boundary conditions 
which we already mentioned in section 4. Of course we have enough integration 
constants to satisfy all the boundary conditions. Especially, for the scalar field, we 
set its value at the horizon to be 0o and at infinity it takes any free value. This 
proves the attractor behavior of the system for scalar field moduli. 

The re-expansion of (p^ (r) shows that it includes terms like x n log x which indi- 
cate the non-analytic behavior, though regular, near the horizon. To avoid quoting 
lengthy results, we demonstrate our results in figure 2. 
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Figure 2: 0(r) vs. log(r), where numerical coefficients for the potentials are #2 = 
Vi>2 = 1 and gi = 2. Different curves represent different asymptotic values for (p^. 
The attractor point is 0o = at the horizon, r# = 1. 

6 Conclusion 

In this paper, we studied the attractor mechanism in a theory of gravity coupled 
to gauge fields and scalar fields, with higher derivative R 2 terms in the action. 
By investigating solutions of the equations of motion, we observed the attractor 
behavior explicitly. We looked for all possible solutions which admit the criteria of 
being regular at the horizon and free in the asymptotic region. 

The near horizon analysis done in section 3 (see also [12]), shows the criteria 
for attractor behavior. Although, the existence of a consistent set of solutions in- 
volving the backreaction remains far from obvious. In R 2 gravity and especially in 
Gauss-Bonnet theory, we solved the equations for the background and scalar fields 
simultaneously which means inclusion of backreaction effects. 

For the Gauss-Bonnet theory, given flat asymptotic boundary condition for back- 
grounds, the regularity of scalar fields at the horizon is a sufficient (and obviously 
necessary) condition to meet the attractor mechanism. The regularity at the hori- 
zon, in the non-analytic sector, sets some restriction on the effective potential Wh 
which is indeed its minimality condition at its critical point. 

The solution with analytic behavior at the horizon appears as a new branch in 
the set of attractor solutions which gives no minimality restriction on the effective 
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potential, Wh- So we observe that, compared to the case studied in [6], turning on 
the R 2 terms and coupling them to the moduli in the form of ()29|). a new sector 
of solutions appear for which the extremality of the effective potential is enough to 
produce the attractor behavior. 

Although, for the analysis in this paper we considered asymptotically flat space 
times, with the same technology, it is easy to check that attractor mechanism works 
as well for asymptotically (A)ds backgrounds. 
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